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Abstract 

In this paper, a study of light-cone distribution amplitudes (LCDAs) for s-wave heavy meson 
are presented in both general and heavy quark frameworks. Within the light-front approach, 
the leading twist light-cone distribution amplitudes, 4>m(u), and their relevant decay constants of 
heavy pseudoscalar and vector mesons, /m, have simple relations. These relations can be further 
simplified when the heavy quark limit is taken into consideration. After fixing the parameters that 
appear in both Gaussian and power-law wave functions, the corresponding decay constants are 
calculated and compared with those of other theoretical approaches. The curves and the first six 
^-moments of 4>m(u) are plotted and estimated. A conclusion is drawn from these results: Even 
though the values of the decay constants of the distinct mesons are almost equal, the curves of their 
LCDAs may have quite large differences, and vice versa. Additionally, in the heavy quark limit, 
the leading twist LCDAs, &Q q (uj) and $>Q q (u), are compared with the 5-meson LCDAs, ip + (uj), 
suggested by the other theoretical groups. 
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I. INTRODUCTION 



The decay constants of heavy mesons with a c or a b quark are significant quantities and 
they play an important role in studies of CP violation, Cabibbo-Kobayashi-Maskawa (CKM) 
matrix elements, the D — DoiB — B mixing process, and leptonic or nonleptonic weak decay. 
Experimentally, new data on the decay constants of the pseudoscalar mesons fo and fo s have 
been reported l|-|4j which has provided a precise method for comparing different theoretical 
calculations and for checking their accuracy. During the last decade, the decay constants 
of both pseudoscalar and vector heavy mesons have been studied by lattice simulations j^J, 
the relativistic quark model js-lO], and the field correlator method 11]. The light-cone 
distribution amplitudes (LCDAs) of hadrons are key ingredients in the description of the 
various exclusive processes of quantum chromodynamics (QCD), and their role is analogous 
to that of parton distributions in inclusive processes. In terms of the Bethe-Salpeter wave 
functions, <p(ui,ki±), the LCDAs, are defined by retaining the momentum fractions, 

Ui, and integrating out the transverse momenta, ki± 12|. They provide essential information 
on the nonperturbative structure of the hadron for the QCD treatment of exclusive reactions 
and they play a central role in all known factorization formulas. Specifically, the leading 
twist LCDAs describe the probability amplitudes for finding the hadron in a Fock state with 
the minimum number of constituents. Experimentally, the fact that I?-physics exclusive 
processes are under investigation in BaBar and Belle experiments also urges the detailed 



study of hadronic LCDAs. In the literature, the LCDAs of heavy quar 
estimated by various nonperturbative approaches, such as QCD sum rules 



ionia have been 
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factorization 



16|, NRQCD 



171 ]. and the light-front quark model 18|-|20|. As for heavy-light mesons, the 



CDAs of .B-meson %p± were first introduced within the heavy quark effective theory (HQET) 



2jJ , and the following studies were intensive 



22 



whereas the ones of other heavy-light 



37J. 



mesons were discussed in a non-HQET framework 

In the past decade, the most significant progress made in the QCD description of hadronic 
physics was, perhaps, in the avenue of heavy quark dynamics. The analysis of heavy hadron 
structures has been tremendously simplified by the heavy quark symmetry (HQS) proposed 
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sgur and Wise [38|, |39[ and HQET developed from QCD in terms of 1/mg expansion 



42]. HQET has provided a systematic framework for studying symmetry breaking I/tuq 



corrections (for a review, see 43[). Moreover, in terms of heavy quark expansion, HQET 
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offered a new framework for the systematic study of the inclusive decays of heavy mesons 



44j-|47l|. However, the general properties of heavy hadrons, namely their decay constants, 



transition form factors and structure functions etc., are still incalculable within QCD, even 
in the infinite quark mass limit with the utilization of HQS and HQET. Hence, although 
HQS and HQET have simplified heavy quark dynamics, a complete first-principles QCD 
description of heavy hadrons has still been lacking due to the unknown nonperturbative 
QCD dynamics. 

This paper has focused on the study of the decay constants and the leading twist LCDAs 
of pseudoscalar (D, D s , B, B s , B c ) and vector (D*, D*, B*, B*, B*) mesons within both 
general and heavy quark frameworks. From the definitions of the decay constant and LCDA 
(or quark distribution amplitude (DA)) |l2|, these two properties seemed to be closely re- 
lated. In terms of a detailed analysis, the purpose of this study is to transparently realize 
the relation between the decay constant and LCDA of the heavy meson. We believe that a 
thorough understanding of these universal nonperturbative objects would be of great benefit 
when analyzing the hard exclusive processes with heavy meson production or annihilation. 
Additionally, in this study, the s-wave heavy meson has been explored within the light-front 
quark model (LFQM), which is a promising analytic method for solving the nonperturbative 
problem of hadron physics [48j] , as well as offering many insights into the internal structures 
of bound states. The basic ingredient in LFQM is the relativistic hadron wave function 
which generalizes distribution amplitudes by including transverse momentum distributions, 
and which contains all the information of a hadron from its constituents. The hadronic 
quantities are represented by the overlap of wave functions and can be derived in principle. 
The light-front wave function is manifestly a Lorentz invariant, expressed in terms of internal 
momentum fraction variables which are independent of the total hadron momentum. More- 
over, the fully relativistic treatment of quark spins and center-of-mass motion can be carried 
out using the so-called Melosh rotation Q|. This treatment has been successfully applied 
to calculate phenomenologically many important meson decay constants and hadronic form 



factors 



50 



55[. Therefore, the main purpose of this study was the calculation of the leading 



twist LCDAs of s-wave heavy mesons within LFQM. 

The remainder of this paper is organized as follows. In Sec. II, the leading twist LCDAs 
of s-wave heavy meson states are derived within general and heavy quark frameworks. In 
Sec. Ill, the formulations of LFQM within the general and heavy quark frameworks are 
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reviewed briefly, and the decay constants and the leading twist LCDAs then extracted. In 
Sec. IV, numerical results and discussions are presented. Finally, the conclusions are given 
in Sec. V. 



II. LEADING TWIST LCDAS OF s-WAVE MESONS 
A. General Framework 

The amplitudes of the hard processes involving s-wave mesons can be described by the 
matrix elements of gauge-invariant nonlocal operators, which are sandwiched between the 
vacuum and the meson states, 

(0\q(x)T[x,-x]q(-x)\H(P,e)), (2.1) 

where P is the meson momentum, e is the polarization vector (of course, e does not exist in 
the case of pseudoscalar meson), T is a generic notation for the Dirac matrix structure, and 
the path-ordered gauge factor is: 



[x, y} = P exp 



ig s / dt(x - y)^(tx + (1 - t)y) 



o 



(2.2) 



This factor is equal to unity in the light-cone gauge which is equivalent to the fixed-point 



gauge, (x — |/) M yl M (x — y) — 0, as the quark- ant iquark pair is at the lightlike separation [56 ]. 
For simplicity, the gauge factor will not be shown below. 

The asymptotic expansion of exclusive amplitudes in powers of large momentum transfer 
is governed by the expanding amplitude, Eq. (12. ip . shown in powers of deviation from the 
light-cone x 2 = 0. The two lightlike vectors, p and z, can be introduced by 

p 2 = 0, z 2 = 0, (2.3) 



so that p — > P in the limit Mfj — > and z x for x 2 = 0. From this it follows that 57] 

1 



M 2 H 



Px - J(Px) 2 - x 2 M 2 H 



= ^-P^ + 0(x'), 
M 2 

= P" - (2.4) 
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where Px = P ■ x and Pz = pz = y (Px) 2 — x 2 Mfj. In addition, if the meson is assumed to 
move in a positive £3 direction, then p + and z~ are the only nonzero components of p and z, 
respectively, in an infinite momentum frame. For the vector meson, the polarization vector 
e M is decomposed into longitudinal and transverse projections as 



— [p^-z* 
pz \ 



M 



H 



2pz 



e" - eff, 



(2.5) 



respectively. 

LCDAs are defined in terms of the matrix element of the nonlocal operator in Eq. ( 12. ip . 
For the pseudoscalar (P) and vector (V) mesons, LCDAs can be defined as 



(0\q(z)Yl 5 q(-z)\P(P)) = ifp / due*" 



P^4>p{u) + 2? 



2pz 



9pW 



(0\q(z)Yq(-z)\V(P,e x=0 )) = f v M v jf du e ^ z {^—<h\\ («) + ^ ± 9v±(' 

;M*g V3 (u)}, 



2(pz) 



(2.6) 



(2.7) 



(0\q(z)a^q(-z)\V(P,e x=±1 )) = ft f du e^{(e^ - e>")0y» 



Mlez . 



+ {j?z v - p v z^X-P^h v \\ (u) + (ef z v - e\ 



v 



{pz) 2 ' _J -" '2/c 

where it is the momentum fraction and £ = (1 — w) — it = 1 — 2m. Here 0p, and 
<pv± are the leading twist-2 LCDAs, and the others contain contributions from higher-twist 
operators. The leading twist LCDAs are normalized as 



du(f>(u) 



and can be expanded 58( in Gegenbauer polynomials C^' 2 (£) as 



E^)Q 3/2 (0 



(2.9) 



(2.10) 



where as (O = 3(1 — £ 2 )/4 is the asymptotic quark distribution amplitude and ai{p) are 
the Gegenbauer moments which describe to what degree the quark distribution amplitude 
deviates from the asymptotic one. Cf (£)' s have the orthogonality integrals 



2(/ + l)(Z + 2) 



Sim- 



-l- 2/ + 3 

Then a; can be obtained by using the above orthogonality integrals as 



(2.11) 



An alternative approach to parameterize the quark distribution amplitude is to calculate 
the so-called ^-moments: 

= f^dZCM,!*), (2.13) 

as calculated in this work. 

To disentangle the twist- 2 LCDAs from higher twists in Eqs. (12. 6 j) ~ (12.8)) . the twist-2 
contribution of the relevant nonlocal operator q(z)Tq(— z) must be derived. For the V = 
7^(75) case, the leading twist-2 contribution contains contributions of the operators which 



are fully symmetric in the Lorentz indices 
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60 : 



1 f(z-D) n niz-D)" 1 ' 1 ^- ~\ 
[q(-z)n l5 )q(z)] 2 = £ -g(0){^— f r + U{Z J £>V (75)g(0), (2.14) 

where D = D — L) and — igB a (\ a /2). The sum can be expressed in terms of a 

nonlocal operator, 

[q(-z)-f( l5 )q(z)] 2 = C dt^-q(-tz) }t(r( 5 )q(tz). (2.15) 

Jo oz^ 

Taking the matrix element between the vacuum and the s-wave meson state, we obtained: 
(0\[q(-z)Yl 5 q(z)h\P(P)) = ifp £ du^p^Lfe*?' + (P» - p>>) jf* rfte^ 2 |(2.16) 
-z)Yq(z)h\V(P,e x=0 )) = f v M v du<p v \\(u)L»yS pz 



+ V ~ ^~z) Jo dte ^ PZ \ ( 2 - 17 ) 



The derivations of (I2.17P as shown in Ref . [59| , are applied to those of Eq. (I2.16P . We use 
Eq. (I2.14p . and then expand the right-hand sides of Eqs. (12. 16)) and (12. 1 T)) as 

°° 1 /n\-/n\{( z -D) n u n(z-D) n ~ l 



E ^<0|9(0){^-y + n{Z '" " D»A l5 q(0)\P(P)) 
^ n\ I n + 1 n + 1 J 

00 -n r r 1 1 

% fpT,- duMmpzrU + (p"- P n ^H, (2.1s) 

n=0 n - JC) I J0 J 

E ^(01^(0)1^^7- + n(g, ^ 1 S>»4 g (o)|y(p >e )) 

„ „n! L n + 1 n + 1 J 



/ V M V £ ^ / 1 ^||(«)(^) n (^- + fe" /Wl, (2.19) 

^ n! 7o [ pz \ pz J Jo J 
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respectively. Picking n = in Eqs. (I2.18P and (I2.19p . we obtain 



(0|g(0) 7 M 7 5 g(0)|P(P)) = if P P» / duMu), 



(0|g(0)7^(0)|^(P),e A=0 ) = f v M v e^ du(f) n (u) 

Jo 

From the normalization of Eq. (12.91) . we have 

(0|g7" T5 g|P(P)> = if P P^ 
(0\q^q\V(P,e)) = fvM v e» 



(2.20) 
(2.21) 



(2.22) 



which are taken as the definitions of decay constants fp and fy in the literature. 

Next, we consider the case of T = a^, where the leading twist-2 contribution contains 
contributions of the operators: 

[q(-z)a^q(z)} 2 = ^ lg(Q) j^L^ + ^ U \ 



n=0 



2n+l 2n + l 

n(z ■ D) n ~ l 



2n + l 

The sum can also be represented in terms of nonlocal operators 



DV" g(0). 



(2.23) 



iq(-z)a^q(z)] 2 = fdt 
Jo 



<} q(~t 2 z)a' u q(t 2 z) + z a ^-q(-t 2 z)a^ a q(t 2 z) 



dz 



dz v 



(2.24) 



Returning to Eq. (12 .8p . it can be rewritten as 

(0|9(*)<V?(-*)I^<*= ±1 )) = & f due^ z {^P v -e v P^ VL (u) 

J 



+ (PpZv-PuZp) ,jll [ h v\\{u) - <j)v±{u)\ 



(pz) 2 
2pz 1 



(2.25) 



We sandwich both sides of Eq. (12.241) between the vacuum and the vector meson state as 

[z)a„ v q(-z)]2\V(P,e x=±1 )) 
^\q(-t 2 z)a. v q{t 2 z)\V{P,e))+z a ^m 

pi 



dt 











<j>v±(u) 


(e^Pu - e v P^) f 






J 



+(h v \\{u) - <fr v ±(u)) 



J dte^ t2pz + 2 V z%M) J q dtt 2 e^' 



pz 



. (2.26) 



The integral is performed as 

- 1 - — l r 1 , d , 



£ / dtt 2 e^ pz = / dtt^-e i&2pz = e^ pz - [ dte^ t2pz 

Jo 2pz Jo at 2pz L Jo 



10 2pz Jo at 2pz 

and then we substitute Eq. (12.261) for Eq. ( 12.271) to obtain 

(0\[q(-z)*^q(z)) 2 \V(P,e x=±1 )} 



[2.27) 







[du< 


(f>v±(u) 







where 



S 



+ (h V \\(u) - <f>v±{u) 



M 2 

y ' v L 1 ; 2pz 



(2.28) 



V 



f7~M 2 

v (v^z v - v v z^ 
2{pzf KV P > 



M 2 

W v = —^(e"z u -€ u z tl ) 
pz 



(2.29) 



In contrast to Eqs. H2. 161) and (12. 17ft . the twist- 2 LCDAs do not disentangle entirely from 
the higher twists in Eq. (12. 28ft . Taking the product with e± fl z u in Eq. ( 12.28ft to obtain 

z)^'e ± ^ 5 q(z)] 2 \V(P,e x=±1 )) = ift [' du0 v± (u)he ■ e ± Pz) \e^ z + P dte****\, 

Jo 2 L Jo 

(2.30) 



we then use Eq. (I2.23P and expand the right-hand side of Eq. (I2.30p as 

E^(0 1 (0) {n \ 1){ *: Dr ^'e^ 5q (0) | V(P, e A=±1 ) > 
ZTk n\ 2n + 1 



n=0 







1 -f 


f 1 dtt 2n 




Jo 



00 i n r 1 1 
= /vEl du<t>v±{u)-{e-e L Pz){&zy 

n=0 U - J ° 1 

Picking n = in Eq. (12. 31ft . we obtain 

(0|g(0K'e ±M g(0)|V(P,e A=±1 )> = # C du<p v± (u)(e ■ e ± Pz). 

Jo 

From the normalization of Eq. ( 12.91) . we have 

(0|g(0)^-e^g(0)|V(P,e A=±1 )) = f^(e-e ± Pz), 
which is consistent with the usual definition of fy as 

(0|g(0)tr^(0)|y(P,e A=±1 )) = ft{?P v - <?P»), 
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(2.31) 



(2.32) 



(2.33) 



(2.34) 



B. Heavy Quark Framework 



In general, the theoretical description of meson properties relies on the bound state models 
with a relativistic normalization: 



(M(P')\M(P)) = 2P (27^) 3 ( 5 3 (P , - P). 



(2.35) 



At low energies, however, these models have little connection to the fundamental theory of 
QCD. Then the reliable predictions are often made based on symmetries. A well-known 



example is HQS 



43], which arises since the Compton wave-length, l/m,Q, of a heavy quark 



bound inside a hadron is much smaller than a typical hadronic distance (about 1 fm), and 
rriQ is unimportant for the low energy properties of the state. For a heavy-light meson 
system, it is more natural to use velocity instead of momentum variables. Then it is 
appropriate to work with a mass independent normalization of a heavy-light meson state: 

(M(v')\M(v)) = 2v°{2n) 3 5 3 {Av' - Av), (2.36) 

where A = M — itlq is the so-called residual center mass of a heavy-light meson. The relation 
between these two bound states is 

(2.37) 



\M{P)) = VM\M(v)). 



In addition, the heavy quark field can be expanded as 

1 



43] 



Q(x) 



-irriQV-x 



1 + 



h v (x), 



(2.38) 



iv ■ D + Ivtiq — is 

where h*(x) is a field describing a heavy antiquark with velocity v. Then the current qTQ 
can be represented as 



h„. 



(2.39) 



Substituting Eqs. (1237]) and fl2T39|) into the definitions of LCDAs, Eqs. dHJ ~ flZU) give 



(0\q(z)YlM-z)\P(v)) = iF P / du e % 

Jo 



2vz 



(0\q(z)rh(-z)\V(v,e x=Q )) = F v du e^ z U— + e£GV±(w) 

Jo 1 vz 



ez 



jGV 3 (w)}, 



(0\q(z)a^h v (-z)\V(v,e x=±1 )) = Fy 



2(vz) 2 

1 r 

du e lu)VZ {(elv v - e u ± v^ v± (oo) 



(2.40) 



(2.41) 



ez I 
+ (v»z v -v v z»)——;H v Auj) + (e»z v -e^z^—H 



(vz 



2vz 



vs 



where Fm = vM/m, = 4>i(u)/M } and u was first introduced in Ref. [6l|] as the 

product of longitudinal momentum fraction u of the light (anti) quark and the mass of heavy 
meson M, namely u = uM . Following a similar process, the leading twist LCDAs are 
obtained as 



(0|g(0) 7 M 75^(0)|P(i;)) = iF P v» / dw$ P (u) 



(0\q(0)^h v (0)\V(v),e x=o ) = F v e» / du$ vll (uj) 

Jo 





oo 



(0\q(0)v' M e-L»hv(0)\V(v,ex=±i)) = F^(e-e ± vz) / du$ v± (tu). (2.45) 

Jo 



v\\ 

oo 



(2.43) 
(2.44) 



The authors of Ref. 



2l| defined two quark-antiquark wave functions in momentum space 



i/j±(lu) of a heavy-light meson in terms of the matrix element: 



(0\q{z)Th v {-z)\M{v)) = f / e^duTr 



where / = -Fm/2 and 



M(v) 



V+(w) + 2^[^-(w)-^+(w)l 



1+ i> [ ~~ Z T5' ^ or pseudoscalar meson M(v) 



7W(t;)r (2.46) 



(2.47) 



2 I for vector meson M*(v,e). 

Evaluating the trace for various choices of T and taking the heavy quark limit, they obtained 



= $yy(w) = $vu.(w) = ^+(w), 



(2.48) 



and the normalization conditions 



dcji[) + (cj) = 1. 



In addition, the authors of Ref. 



211 ] defined the moments of ip+(uj) as 



(2.49) 



(2.50) 



and used the equations of light and heavy quarks to obtain the relation between the first 
moment and the residual center mass: 



(">+ = |a. 



(2.51) 
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III. FORMULISM IN LIGHT-FRONT APPROACH 



A. General Framework 

An s-wave meson bound state, consisting of a quark, q 1 , and an antiquark, q 2 , with total 



momentum P and spin J, can be written as (see, for example 



52|) 



\M(P, S, S,)) = J {d 3 h}{d 3 k 2 } 2(2vr) 3 5 3 (P - h - k 2 ) 



X £ V ss *(k u k 2 ,\ u \ 2 ) \q l (k u X 1 )q 2 (k 2 ,X 2 )}, (3.1) 



where k± and k 2 are the on-mass-shell light-front momenta, 



and 



k = (k + ,k ± ), k ± = (k\k 2 ), k-= mq k+ ± , (3.2) 



{d A k} \ 



2(2tt) 3 ' 

|g(Ari,Ai)g(fc2,A 2 )) =6t(fci,Ai)dt(A; 2s A 3 )|0), (3.3) 
{b(k r , A'), 6 f (A;, A)} = {d(k', A'), rf f (fc, A)} = 2(2tt) 3 5 3 {k' - k) 5 X >\. 

In terms of the light-front relative momentum variables (u,k±) defined by 

k+ = (1 - U )P+, fc+ = «P + , 

fell. = (1 - u)P± + K±, k 21 _ = uP± - K_l, (3.4) 

the momentum-space wave function ty SSz can be expressed as 

y ss *(h, h, Ai, A 2 ) = -L^-RSS-fa K± ) cpfa K± ), (3.5) 

viv c 

where <p(u, k±) describes the momentum distribution of the constituent quarks in the bound 
state, and Rf^L constructs a state of definite spin (S, S z ) out of the light-front helicity 
(Ai,A 2 ) eigenstates. Explicitly, 

R M\ 2 ( U , K ±.) = Yl ( a iI 7 ^a/( 1 - M > Kx,m 1 )\s 1 )(X 2 \TZi I (u, -K±,m 2 )|s 2 )(- -; sis 2 |--; SS£$.6) 



where |sj) are the usual Pauli spinors and TZm is the Melosh transformation operator 50]: 



(s|/CAr(u»,«±,m»)|A) = — , (3.7) 

(rm + MiM ) 2 + k\ 
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with iii = 1 — u, u 2 = u, and n = (0, 0, 1) a unit vector in the indirection. In addition, 



Ml 



(ei + e 2 ) 2 = -± ^ + - 



1-u 



u 



(3.8) 



mf + k 2 , + nl. 



where k z is the relative momentum in z direction and can be written as 

uMq ml + k\ 



(3.9) 



2 2uM 

M is the invariant mass of qq and generally different from mass M of the meson which 
satisfies M 2 = P 2 . This is due to the fact that the meson, quark and antiquark cannot be 
simultaneously on-shell. We normalize the meson state as 



<M(P', S', S' Z )\M(P, S, S g )) = 2(2vr) 3 P+5 3 (P' - P)6 S , S 6 S , S , 



(3.10) 



in order that 



du d 2 n±_ 



\ip(u,K ± )\ 2 = 1. 



2(2vr) 3 

In general, for any function <p(u, k±) has the form of 



(3.11) 



where normalization factor N is determined from Eq. ( 13. lip . 



(3.12) 



In practice, it is more convenient to use the covariant form of PaiA 2 [50, [54], [62|: 



1 ^2 



y/2 M (M + mi +m 2 



-u(k u X 1 )(P + M )Tv(k 2 ,X 2 ) 



(3.13) 



where 



M = ^M 2 - ( mi - m 2 ) 2 , 
2m 
~k+ 
2m 



u(k, X)u(k, A') = ~^~S\ x 
v(k, X)v(k, A') = ~r~i~$\,\ 
For the pseudoscalar and vector mesons, we have: 



P = k 1 + k 2 , 
u(k, X)u(k, A) 



5>(MMM) 

A 



ft + m 

ft — m 
= k+ ' 



(3.14) 



r P = 75, rv = - /(A), 



(3.15) 
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where 



— r ± (±l) • P ± , 0, r ± (±i) 



t A=±l 

eL(±l) = T(l,±i)/V2, 
M 



P + ,Pi 



(3.16) 



Equations (13.131) and ( 13. 15ft can be further reduced by the applications of equations of motion 



on the spinors 



where 



54|: 



R ss 



Tp = 75, r' v 



u(h, Xi)T'v(k 2 , A 2 ) 



e-(*i 



(3.17) 



(3.18) 



M + mi + m 2 

Next, the matrix elements of Eqs. (J22DJ), (E2H), and (12T321) are calculated within LFQM, 
and the relevant leading twist LCDAs are extracted. For the pseudoscalar meson state, we 
substitute Eqs. (ED), (ETSjl . and (I3TT7D into Eq. (I2T20D to obtain 

(Ol^^sgilW) = N c J{d 3 h} £ * s ^(fc l5 A 1; A 2 )<0lg 2T ^7 5 aiki? 2 > 



Ai,A2 



'k k 

iJN c I {d 3 kA y —J-J ipTr 
z/pP" / du^(u). 



7^1 



kf 



75 | 



For the "good" component, /i = +, the leading twist LCDA, <pp, is extracted as 

2y/E f d 2 K L [(1 — u)m2 + KroJ 



/p J 2(2vr) 



u(l - m)M 



-(f(u, K±). 



(3.19) 



(3.20) 



A similar process is used for the vector meson which corresponded to Eqs. ( 12.211) and ( 12.321) . 
and then the leading twist LCDAs are extracted as 



<f>v\\{u) 



4>vl{u) 



2\/E f d 2 n± <p(u,k±) 



fv J 2(2tt) 3 ^ u (i _ u ) Mq 

2a/6 r d 2 K± (p(u,Kx) 
If' 2(2tt)3 y^TT^Mo 



um\ + (1 — u)m 2 + 
umi + (1 — u)m2 + 



2k\ 



Mq + mi + m 2 

2 



M + m 1 + m 2 



(3.21) 



(3.22) 



From the normalization of Eq. (12.91) . we found not only that the equations of fp and fv 
were consistent with that of [6j, but also that the decay constants and the leading twist 
LCDAs has the simple relations 

f P + fv = 2f$, <j> P (u) + <h\\(u) = 2<j>v ± (u). (3.23) 
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B. Heavy Quark Framework 



If one takes mi = wlq — > oo, that is, the heavy quark limit in the heavy-light meson, 
then two inequalities, rriQ ~ M ^> m 2 , k± and w — > 0, are obtained. From Eqs. (13.201) ~ 
f)3.22p . the decay constants and the leading twist LCDAs are simplified as 

fp - fv - fv oc F M} (f) P ~ fcw ~ (f) V1 _ oc $ M , (3.24) 

which are independent of the form of These are consistent with HQS between 

the s-wave heavy-light mesons. The exact form of <3>m 5 however, must be derived by the 
redefinition of the meson bound state. Let us consider the bound states of heavy mesons in 
the heavy quark limit: 



\M(v;S,S z )) = J{d 3 q}{d 3 k 2 }2(2rt) 3 5 3 (Av-q-~k 2 



x J £^ ss '{u, K± ,X 1 ,X 2 )bl{q,X 1 )S(k 2 ,X 2 )\Q), (3.25) 

Ai,A2 

where q = k\ — vtiqV is the residual momentum of heavy quark. The operators bl(q, X\) 
create a heavy quark with 

{&„(<?, Ai), ft^^Ai)} = 2(27rfS vvl 5 3 (q-q')S XlK , (3.26) 

The relative transverse and longitudinal momenta, k± and obtained by 

u m 2 2 + k 2 ± . , 

k x = k 2± - uv ± , k z = — . (3.27) 

2 ZOJ 

The momentum-space wave-function ty SSz can be expressed as 

q ss *(cu, «l, Ai, A 2 ) = 4r%( w . «J-) ^(^' ^J-) j (3-28) 



where 

^+ 

i? 5& (w, k_l, Ai, A 2 ) = - 2 — u(v, X 1 )Tv(k 2 , A 2 ) (3.29) 

v^V^ + m 2) 2 + «i 

withr = 75 (- ft) for 5 = 0(1), 



e A=±l 



^e ± (±l)-v ± , 0,el(±l) 
3U> = f zl 7^ ; A^, (3.30) 
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and u(v, Ai) is the spinor for the heavy quark, 



J2u(v,X)u(v,X) = ^. (3.31) 

A V 

The normalization of the heavy meson bound states can then be given by 

(M(v', S', S' z )\M{v, S, £7 Z )> = 2(2n fv + 5 3 (Av' - Av)5 S s'$s z s' z , (3.32) 

which not only leads to Eq. (12. 3 7ft . but also to the space part (p SS!1 (uj, k\) (called the 
light-front wave function) in Eq. (13.25]) which has the following wave function normalization 
condition: 

In principle, the heavy quark dynamics are completely described by HQET, which is given 

by the 1/tuq expansion of the heavy quark QCD Lagrangian: 

oo / i \ n 

L = Q{ijp- m Q )Q = ( A, (3.34) 

are then determined by the leading Lagrangian 



n=0 



63| have shown, from the light-front bound state 



Therefore, \M(v;S,S z )) and ^^(u^i) 
Cq = h v iv ■ Dh v . The authors of Ref. 
equation, that (p SSz (U, k,\) must be degenerate for S = and £7=1. As a result, we can 
simply write 

SSz {lo,k 2 ± ) = 0(u,k 2 ± ) (3.35) 

in the heavy quark limit. Equation (I3.25P together with Eqs. f !3.29|) and (I3.35P are then the 
heavy meson light-front bound states in the heavy quark limit that obeyed HQS. From the 
normalization conditions of Eqs. (13. lip and (13.331) . we obtain the relation between the wave 
functions ip(u, n 2 ± ) and <p(oj, k^_): 

cp(u, k 2 ± ) = VM0(cu, k 2 ± ). (3.36) 

Next, the matrix elements of Eqs. (12.431) . (I2.44p . and (12.451) can be calculated, and the 
relevant leading twist LCDAs extracted. For the pseudoscalar meson state, we substitute 
Eqs. (153511 . (|53B)» . and <^M> into Eq. <^M) to obtain 

(0|g 7 M 75^|^)> = N c J{d 3 k 2 } £ * SS '(u,K±,Xi,X2)(0\ta7*7sK\Q<h) 



Ai,A2 

J3, 



i ViV c / {d 6 k 2 } fiTr 
V2 J(u + m 2 ) 2 + k\ 



iFpv^ J du$ P (uj). (3.37) 
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For the "+" component, the leading twist LCDA $p is extracted as 



$p(u;) 



2x/6 r d 2 K± 



--{p(uj, K±). 



(3.38) 



Fp J 2(2tt) 3 v /( w + m2 )2 + K 2 ± 

In contrast with (f)(u), $(w) represents the distribution of the longitudinal momentum car- 
ried by the light degree of freedom. A similar process is used for the vector meson which 
corresponds to Eqs. (12.441) and (I2.45H . and the results are 



F P = F v = F$, $p(w) = = $y±(w), 



(3.39) 



which are consistent with Eqs. (I2.48j) and (I2.49p . 



IV. NUMERICAL RESULTS AND DISCUSSIONS 

In this section, the decay constants and LCDAs of D^*\ D^*\ B^*\ Bf\ and B^ are 
studied. We consider two kinds of J-'(IkI), one is the Gaussian type, the other is the power- 
law type: 



J™(\k\) = 
F>{\k\) = 

then the corresponding wave functions are 



exp 



2f3 2 , 
1 ^ 



(it, Kj_ J 



LfP{U,K ± ) = a I — 



eie 2 



4(^ 3/4 / 



\k\ 2 /^ 



exp 



(4.1) 
(4.2) 



o 



2mM / 



2/3 2 



27r\V 2 / 



eie 2 



w)M 



£ 5 



1 2 



Kl + 



(4.3) 
(4.4) 



2 2«M 

and can be used to calculate decay constant /, the LCDAs (f)(u), and the ^-moments of 
<p(u). Prior to the numerical calculations, the parameters mi, m 2 and /3, which appeared 
in the wave function, have to first determined. For the light quark masses, we used the 
decay constants /„-, fx and the mean square radii (r 2 +), {r 2 K o) to fit m Uj d{= m q ) and m s 



64|. For the heavy quark masses, however, the relevant measurements are insufficient. We 



determined m c and m;, by the mass of the spin-weighted average of the heavy quarkonium 
states and its variational principle for the relevant Hamiltonian {20 1. 

As regards parameter (3, it is determined by the decay constant of the heavy meson. 
Recently the CLEO collaboration updated their data concerning Br(.D + — > and an 
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average value was reported |65j: fo+ = 206.0 ± 8.9 MeV. In addition, the authors of Ref. 



66| averaged Br(B — > r v) from the Belle |67j and Babar 68|, [69] collaborations and 
extracted fs = 204 ±31 MeV. The parameters (3 cq and f3b q can then be determined. As 



mentioned in the previous work 
SU(3) symmetry breaking, that is, m s /m q as follows: 

3 



64j, the ratios, (3 CS / (3 cq and Pbs/ Pbq can be related to the 



AM 



D a D* 



AM, 



B 3 Bi 



'A 



AMr 



(4.5) 



1 dd* rn s \(3cqJ AM BB * m s \(3 bq J 

Therefore, (3 CS and fibs are not independent parameters. Concerning the decay constants 
of B c , we quote the average result of QCD sum rules [to | : fs c = 360 MeV to extract the 
parameter /?&,. All the parameters are listed in Table 1. 

TABLE I: Input values of quark masses and /3's (MeV). 





m q 


m s 


m c 


m b 




Pes 






Pbc 




251 


445 ± 36 


1380 


4780 


465 ± 22 


567 ± 42 


587 ± 74 


727 ± 120 


815 




172 


296 ± 12 


1360 


4770 


505 ± 25 


608 ± 40 


575 ± 77 


706 ± 113 


815 



Next, we used the parameters in Table 1 as input to calculate the decay constants fp, 
fv, and fy of the relevant heavy mesons. The values of the ratios fv/ fp, fp' / fp, fv / fv 

are also included. Tables 2 and 3 show a comparison of the results of this work with other 

I I 

theoretical calculations. In a previous work [64j, we pointed out that ratios foj fo and 
fp-sl 1 fp were not only chiefly determined by the ratio of light quark masses, m s /m q , or 
the SU(3) symmetry breaking, but also insensitive to the heavy quark masses m c ^- This 
phenomenon also appears in the ratios /d*//d* and /b*//b* here for both Gaussian and 
power-law wave functions. On the contrary, as shown in Table 3, the ratio fv/ fp is not only 
dependent on the heavy quark mass, but also insensitive to the light quark mass. The reason 
is that, making a comparison between Eqs. ( 13 .20 j) and (I3.2ip . the difference between f P and 
fv is proportional to 2k\/{Mq + m\ +m2). In the Gaussian (power-law) wave function, the 
mean square value of the transverse momentum is equal (proportional) to the square value of 
the parameter (3, or (k^) = (3 2 ((k 2 l ) oc (3 2 ), so the ratio fv/ fp is influenced by the parameter 
(3 and the quark mass. In the case of the different heavy quark, for example, /d*//d and 
/b*//b, as mj is much greater than m c , this effect is greater than that of (3b q > (3 cq , so 
is smaller than On the other hand, in the case of the different light quark, 
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TABLE II: Decay constants of the pseudoscalar and vector heavy mesons (MeV). Linear and HO 
are the different potentials in Refs. [f], 3], FC is the field correlators, BS is the Bethe-Salpeter 
equation, and RQM is the relativistic quark model. 





Experiment 


This work a 


Linear(HO) 


fc[ii] 


BS[8, 9] 


Lattice [51 


RQM [10] 


fn 


206.0 ± 8.9 [65] 


206.0 ± 8.9(206.0 ±8.9) 


211(194) 


210 ± 10 


230 ± 25 


211 ± 14 + ? 


234 


fn* 




259.6 ± 14.6(306. 3t\r$) 


254(228) 


273 ± 13 


340 ± 23 


245 ± 20t ° 


310 


•1 L) 




232.7 ± 11. 7(256. 2+jH) 












fD 

J S 


260.7 ± 6.5[65] 


267.4 ± 17.9(259.7 ± 13.7) 


248(233) 


260 ± 10 


248 ± 27 


231 ± 12±g 


268 


fo* 




338.7 ± 29.7(391.0 ± 28.9) 


290(268) 


307 ± 18 


375 ± 24 


272 ± 12 + 2 n 

— 20 


315 


fn* 




303.1 ± 23.8 (325.3 ±21.5) 












Ib 


204 ± 31 6 


204 ±31(204 ±31) 


189(180) 


182 ± 8 


196 ± 29 


179 ± 18tg 6 


189 


fB* 




225 ± 38(249^2) 


204(193) 


200 ± 10 


238 ± 18 


196 ± 24^^ 


219 


J B' 




214 ± 34(226 ± 37) 












fB. 




281 ± 54(270 ±47) 


234(237) 


216 ±8 


216 ±32 


204± 16±o 8 


218 






313 ±67(335 ±68) 


250(254) 


230 ± 12 


272 ± 20 


229 ± 20±?6 


251 


fk 




297 ±61(302 ±58) 












f Be 




360(360) 


377(508) 


438 ± 10 


322 ± 42 






f B: 




387(423) 


398(551) 


453 ± 20 


418 ±24 










374(392) 













°The value is obtained by ip 9 (ip p ). 

6 This value is extracted by the branching ratio: B(B~ 



t-v) = (1.42 ±0.43) x 10" 



66] 



for example, /d*//d and /dj/Zd,, as m s is slightly greater than m q , this effect is less than 
that of P cs > P cq , so f D */ ' f Ds is a little larger than /d*//d- 

The quark distributions of the heavy meson, (f)p(u), 4>v\\(u), and aV_i_(w) are plotted in 
Fig. 1 and 2. Clearly the difference in the constituent quark masses is greater, the location 
where u peaked is closer to zero. This indicates, relatively, that the lighter the quark, the 
smaller its momentum fraction. We also find that, even though the difference between /d 
and fr>* was more than 25% (almost 50%) for the Gaussian (power-law) wave function, all 
curvilinear distinctions between (pv\\{u) and c J )p(w) are quite small. The reason is that, after 
the k± integration, the curve of 4>m {u) is influenced only by the quark mass, parameter /3, 
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TABLE III: Ratio of the decay constants. In this work, fy/fp = (1 + fy/fp)/2. 





IjA LJl^l 111H_.11 b 


±. 111S WU1 ft. 


_I_jlllt,£lI 1 11W 1 




RSfS Ql 


J_Jd Li Lillet, (J 


ROMh 01 


Id'/Id 




1.26 ±0.02 


1.20(1.18) 


1.27 ±0.05 


1.48 ±0.26 




1.32 


Id* 1 Jd 3 




1.27 ±0.03 
ii.ji m u.uo j 


1.17(1.15) 


1.17 ±0.04 


1.51 ±0.26 




1.18 


Jd 3 1 Id 


1.27 ±0.06 [65] 


1.30 ±0.04 


1.18(1.20) 


1.24 ±0.03 


1.08 ±0.01 


1.10 ±0.02 


1.15 


fn*Jf D * 




1.30 ±0.05 
yi.ZiO m u.uo j 


1.14(1.18) 




1.10 ±0.06 


1.11 ±0.03 


1.02 


fB'/ Ib 




1.10 ±0.02 
^i.zz m u.uo j 


1.08(1.07) 


1.08 ±0.04 


1.21 ±0.27 




1.16 


/b« / /b s 




1.11 ±0.03 

^i.Z^i in U.UO J 


1.07(1.07) 


1.07 ±0.04 


1.26 ±0.28 




1.15 


JB S 1 JB 




l.oo ± U.U i 
(1.32 ±0.08) 


l.Z4(^l .oz J 


1 1 O -U n HQ 

l.iy ± U.Uo 


1 in i c\ ni 
1 .1U ± U.U1 


1.14 zb U.Uo_q 01 


1 1 K 


fBt/fB* 




1.39 ±0.08 
(1.35 ±0.08) 


1.23(1.32) 




1.14 ±0.08 


1.17±0.04t^ 


1.15 


!b* 1 /s c 




1.08(1.18) 


1.06(1.08) 


1.03 ±0.03 


1.30 ±0.24 







°The value is obtained by ip 9 (ip p ). 



and the total spin (that is, the pseudoscalar or the vector meson). As the quark mass and /3 
are the same in (j)v\\(i) {u) and (pp(u), the distinctions between them were slight. On the other 
hand, even though f D * is almost equal to fu a , as shown in Table 2 (or fo'/fo — IdJ Id) 
for the Gaussian wave function, the curvilinear distinction between <j) 9 D *^(u) and <P 9 d{u) is 
obviously smaller than that between 4> 9 d b {u) and 4> 9 D (u). As for the power-law wave function, 
the situation is inverse. Therefore, we can infer that, even though the values of /m's are 
almost the same between the distinct heavy mesons, the curves of 4>m{u) may have the quite 
large differences, and vice versa. Finally, the quark distribution function is displayed in 
terms of the ^-moments, as in Eq. (12.131) . The first six ^-moments (n > 0) are listed in 
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0.0 0.2 0.4 0.6 0.8 1.0 

u 




FIG. 2: Quark distribution amplitudes of the heavy meson for a power-law wave function. The 
solid, dotted, and dashed lines correspond to 4> p p {u), c/jyJu), (j)y ± (u), respectively. 



Table 4. 

For the heavy quark framework, some models for B meson LCDAs have also been adopted 
in the literature. Inspired by the QCD sum rule analysis, the authors of Ref. 
a simple model: 



21] proposed 



(4.6) 
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TABLE IV: First six ^-moments of the s-wave heavy meson. 





(?) 


(?) 


(?) 


(?) 


(?) 


(?) 


rD 


-0.288(- 


■0.251) 


0.210(0.235) 


-0.125(- 


-0.115) 


0.0960(0.111) 


-0.0695(- 


-0.0673) 


0.0558(0.0664) 




-0.213(- 


■0.207) 


0.183(0.217) 


-0.0890(- 


-0.0905) 


0.0738(0.0970) 


-0.0468(- 


-0.0507) 


0.0388(0.0550) 




-0.617(- 


■0.531) 


0.425(0.398) 


-0.312(- 


-0.288) 


0.240(0.234) 


-0.191(- 


■0.185) 


0.156(0.157) 




-0.549(- 


-0.486) 


0.359(0.359) 


-0.254(- 


-0.249) 


0.189(0.200) 


-0.147(- 


-0.154) 


0.117(0.129) 


4: ] 


-0.536(- 


-0.368) 


0.227(0.230) 


-0.133(- 


-0.123) 


0.108(0.0867) 


-0.0553(- 


-0.0527) 


0.0378(0.0403) 



Additionally, the authors of Ref. [25[ suggested a Gaussian-type model: 

,2 



(4.7) 



By applying Eqs. (I2.50p and (I2.5ip . the relation between the residual center mass and the 
parameter A is: 

3 . 



A„ 



2tt 



ii ■ 



(4.8) 



In Ref. |23[, the value Xj = 0.3 GeV corresponded to A q = 0.45 GeV. For a convenient 
comparison, we used this A q and A s = A q + m s — m q to fix parameters f3Q q and Pq s in this 
work. Moreover, the Gaussian wave function ip is given by taking the heavy quark limit in 
Eq. ( 14. 3 p and using the relation Eq. ( 13.36)) : 

7T \ 3 / 4 



<p(u,n±) = 4 



ft 



'1 



m\ + k 2 ± 



exp — 



2w I 



(4.9) 



2u 2 "V 2/3 2 
The light quark masses m q i s \ = 0.251(0.445) GeV are as in Table 1, and we can then obtain 
the values f3n q = 0.279 GeV and (3q s = 0.338 GeV. In terms of these parameters, the leading 
twist LCDAs $>Q q (cu), $q s (o;), ip + i(cu), and ip + u(u) are calculated and plotted as in Fig. 3. 
We find that the curve of &Q q (uj) is close to that of ip + j(u). 



V. CONCLUSIONS 



This study has discussed the leading twist LCDAs of the s-wave heavy meson within the 
light-front approach in both general and heavy quark frameworks. These LCDAs are shown 
in terms of light-front variables and relevant decay constants. In the general frameworks, we 
find that the decay constants and LCDAs of the pseudoscalar and vector mesons have the 
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1.5 



^-N 1.0 

3, 

e 

0.5 
0.0 

1 2 3 4 5 

cd(GqV) 

FIG. 3: Leading twist LCDAs in the heavy quark framework. 

following relations: fp + fv — 2fv an d 0p + 0v|| — 20y±- The parameters m and /3, which 
appear in both Gaussian and power-law wave functions, were determined as follows: (1) the 
light quark masses are fitted by the decay constants and the mean square radii of the light 
meson; (2) the heavy quark masses are determined by the mass of the spin-weighted average 
of the heavy quarkonium states and its variational principle for the relevant Hamiltonian; and 
(3) the hadronic parameter /3's are evaluated by the decay constants of D + , B~, and B c , with 
the former two and the latter one from the experimental data and the average result of QCD 
sum rules, respectively. We find that, for both Gaussian and power-law wave functions, the 
ratios fp*/ fp* and fp*/fp*, as well as fpj fp and fpj 'fp in the previous work, are chiefly 
determined by the ratio of light quark masses m s /m q , or the SU(3) symmetry breaking. On 
the other hand, by making a comparison between fp*/fp, fp*/fp s , fp*/fp, and fp*/fp a , 
the ratio fv/ fp is not only dependent on the heavy quark mass, but also insensitive to the 
light quark mass. 

As shown in Fig.l and 2, we find that even though the difference between fp and fp,* is 
more than 25% (almost 50%) for the Gaussian (power-law) wave functions, all curvilinear 
distinctions between <fiy\\(u) and 4>p(u) are quite small because their main difference come 
from the variations of the quark mass and (3. On the contrary, even though fp* / fp is almost 
equal to fpjfp for the Gaussian wave function, the curvilinear distinction between (j) 9 D ,^(u) 
and <f) 9 D (u) is obviously smaller than that between <f) g Ds (u) and 4> 9 D (u). As for the power-law 
wave function, the situation is inverse. Therefore, we conclude that even though the values 
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of /m's are almost equal among the distinct mesons, the curves of 4>m{u) may have quite 
large differences, and vice versa. 

When the heavy quark framework is used, the above relations for the decay constant and 
the LCDAs can be further simplified as Fp = Fy = Fy and $p = $yn = $yj_, as consistent 
with HQS. For a convenient comparison, the value A q = 0.45 GeV, as suggested in Ref. 23], 
is used to fix /3q 9 (q s ) and to plot the curves of <&q q , $q s , if)+i, and if)+n in Fig. 3. We find 
that the curvilinear distinction between $q 9 (o;) and i/j + i(u) is relatively small. 
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